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1. 


SUMMARY 


A  report,.  PIBAL  number  167,  was  written  by  Salerno  and 

Levine  on  the  buckling  of  cylinders  and  titled  ,fBuckling  of 

Circular  Cylindrical  Shells  with  Evenly  Spaced,  Equal  Strength 

Circular  Ring  Frames.  Part  I.n 

In  discussii  that  report  Dr.  Keunard  pointed  out  certain 

apparent  inconsistencies  in  theory. 

-*■ — ^ 

In  this  present- »iu pert  ^he  material  presented  in  PIBAL 


167  is  discussed  together  with  the  remarks  made  by  Dr.- 


Kenhard,  It  appears  that  some  of  the  theory  presented  in  PIBAL 

f/o, 

numj^Qy- 167  is  open  to  question  but,  nevertheless,  the  conclusions 
obtained  therein  hold  without  modification,  - 


This  report  is  divided  Into  three  main  parts,.  In  part  1 
is  presented  some  material  on  the  theory  ox  instability,  and 
some  of  the  quantities  contributing  to  the  total  potential  energy 
for  the  current  shell  problem  are  derived  in  detail..  In  part  2 
is  discussed  the  essential  differences  between  the  work  of 
Salerno  and  Levine  and  of  Kennard.  It  is  shown  that  although 
some  of  the  work  of  the  former  authors  is  open  to  question  and 
to  criticism  the  terms  representing  the  difference  between  their 
work  and  that  of  Kennard  make  no  essential  difference  to  the  con¬ 
cluding  eigenvalvy formula,  In  part  3  the  effects  of  the  terms 
in  question  in  both  treatments  are  traced  and  considered  in  de¬ 
tail.  As  one  conclusion  it  appears  that  certain  terms  proposed 
by  Kennard,  although  aamis3ible,  lead  to  a  result  different  from 
that  obtained  by  Bryan  for  the  limiting  case  of  an  infinite 
cylinder. . 


LIST  OF  SYMBOLS 

R  o  radius  of  middle  surface  of  shell 

h  »  thickness  of  shell 

L  a  diatajice  between  ring  frames 

L  «  L/R 

o 

V  »  volume  occupied  by  shell 

E  «<  modulus  of  elasticity 

v  «  Poisson’s  ratio 

x,s  **  axial  and  circumferential  coordinates 

9  «  angular  coordinate  circumferentially 

l  »  x/R 

u,v,w  m  displacements  of  shell 

h  **  curvature 

X  **  change  in  curvature 

<7  «■  stress 

e  <=  strain 

"  stress  resultants  per  unit  length 
Mxx*Mss»^xs  c  ,noment  resultants  per  unit  length 
p  *  radial  pressure 

P  »  axial  pressure 

Ua  *  strain  energy  arising  from  stretching  of  middle  surface 

Ub  »  strain  energy  arising  from  membrane  shear  stress 

U  «*  strain  energy  arising  from  bending  moments  and  torques 

induced  in  change  of  state 

Ud  =  strain  energy  arising  from  work  done  by  external  force 
system 


strain  energy  stored  in  the  reinforcing  rings 
work  done  by  radial  pressure 
potential  due  to  axial  pressure 
number  of  half  waves  circumferentially 
number  of  half  waves  axially 
nlB/L 
Kq  «*  ring  parameters 
pR/K 
Ph/K 
h2/12R2 
Eh/(l-v2) 

■=  coefficients  arising  in  expansion  of  determinants 


PART  1.  GENERAL  THEORY 


1.  Instability  Conditions. 

We  consider  a  closed,  right  circular  cylinder  in  an  initial 
state  of  zero  stress.  The  cylinder  is  assumed  bo  be  very  long 
in  comparison  with  its  diameter,  and  to  contain  many  similar  bays 
formed  by  a  set  of  equi-spaced  internal  circular  reinforcing 
rings. 

Let  this  shell  be  now  loaded  by  an  external  hydrostatic 
force  system.  We  desire  to  compute  the  magnitude  of  the  pres¬ 
sure  at  which  the  cylinder  will  buckle  between  adjacent  ring 
frames. 

(1)  Let  the  hydrostatic  load  system  be  applied  to  the  unstressed 
shell  in  such  a  manner  that,  although  on  the  point  of  buckling, 
the  cylinder  has  not  actually  buckled.  Let  the  associated  hydro¬ 
static  pressure  be  p.  .  • 

We  will  call  this  stressed  state  the  state  k .  It  is  an 
equilibrium  state. 

In  the  usual  theoretical  approach  to  a  shell  problem  of  the 
Jr  sent  type  it  is  generally  assumed  that  the  cylinder  is  com¬ 
pressed  uniformly  by  the  hydrostatic  pressure  before  the  rein¬ 
forcing  rings  are  positioned.  Thus  in  state  A  cylinder  generators 
are  straight. 

(ii)  Now,  without  alteration  in  the  magnitude  of  the  hydrostatic 
load  system,  we  enforce  a  small  change  in  the  shell  shape  so  that 
it  takes  on  an  adjacent  bent  shape.  We  will  call  the  associated 
stress-state  state  B.  In  state  B  the  generators  are  no  longer 
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5. 

straignt. 

Such  an  adjacent  deformation  could  be  induced,  for  instance, 
by  the  action  of  a  small  additional  force  system  acting  normal 
to  the  circumferential  wall  of  the  shell. 

If,  cn  the  removal  of  this  small  additional  force  system, 
the  shell  remains  deformed  in  the  adjacent  configuration,  then 
we  say  that  the  hydrostatic  pressure  p  is  the  buckling  pressure 
for  this  shell  structure. 

For  the  buckling  pressure  pf  state  B  as  just  described 
(with  the  small  additional  force  system  removed)  is  also  an 
equilibrium  state. 

(iii)  let  the  total  potential  energy  (i,e.  strain  energy  minus 

work  done  by  external  loads  which  are  assumed  as  constant  in 

going  from  initial  to  final  deflected  shape)  associated  with 

state  A  be  n  . 

a 

Let  the  change  in  potential  energy  in  going  from  state  A 
to  state  B  be  II  . 

D 

Let  the  total  potential  energy  in  state  B  be  n,  measured 
from  aero  stress  content. 

Then,  since  state  B  is  an  equilibrium  state,  the  Potential 
Energy  Theorem  of  Elasticity  states  that 

ft  *  nA  +  iig 

is  stationary.  Hence,  for  small  variations  In  displacements 


oH  «  6(liA  +  %)  -  0. 
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6. 


But  nA  is  also  an  equilibrium  state,  so  that 

611^  «=  0. 

Thus,  for  p  to  be  the  buckling  pressure,  we  require 

6Eb  -  0.  (1) 

(iv)  We  can  obtain  another  result  from  a  slightly  different 
viewpoint.  We  know  state  A  to  be  an  equilibrium  state  and,  for 
the  p  associated  with  state  A  to  be  an  instability  pressure,  we 
^sk  for  the  existence  of  a  neighboring  statu  which  is  also  an 
equilibrium  ?tate  under  the  action  of  the  same  load  system.  For 
.  i  to  be  tne  change  in  potential  energy  in  going  from  state 
A  ro  state  S  must  be  zero.  Hence 

I!b  “  0.  (2) 

We  snail  use  condition  (1)  as  the  criterion  by  means  of 
which  we  will  evaluate,  approximately,  the  buckling  pressure  p. 

We  could  equally  well,  of  course,  use  the  criterion 

6H  *  0. 

2.  Potential  Energy  31. . 

It  is  convenient  to  consider  the  strain  energy  of  t-b=  stress¬ 
ed  shell  as  consisting  of  various  contributions:  those  due  to 
middle  surface  stretching  and  shearing,  and  those  due  to  bending 
and  twisting.  To  do  this  it  is  necessary  to  be  able  to  write 
down  quantities  representing  changes  in  curvature,  etc.,  in  going 
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from  state  a  to  state  B. 

It  would,  perhaps,  be  better  to  consider  the  strain  energy 
in  a  purely  formal  manner,  and  write  it  as 


Be  “  5  ”vaij°iJdT  -  U..1.  -  1.2.3) 


where  the  usual  summation  convention  is  implied,  and 


aij  "  (ei3’ 

via  the  stress-strain  relations.  To  do  this,  however,  would 
require  a  detailed  knowledge  of  all  the  strain  components 
throughout  the  she'll  wall  expressed  in  terms  of  displacement 
components.  This  approach  has  not  been  used  in  PIBAL  number  167, 
and  so  will  not  be  used  here. 

The  potential  energy,  then,  will  be  taken  to  consist  of 
five  contributions: 

(a)  That  arising  from  stretching  of  the  middle  surface  in  going 
from  the  "plane"  state  A  to  the  adjacent  "bent"  state  B,  Stipu¬ 
lating  state  B  to  be  an  adjacent  state  implies  that  it  is  very 
close  to  starve  A  even  though  the  deformation  pattern  is  different. 
As  such  it  is  assumed  that  the  changes  in  magnitude  of  the  normal 
membrane  stresses  are  small  in  comparison  with  their  magnitudes 
in  state  A.  In  fact,  frequently  it  is  assumed  that  the  normal 
membrane  stresses  remain  constant  in  going  from  state  A  to  state 
B.  We  do  not  do  so  here,  however. 

(b)  That  arising  from  membrane  shear  stresses.  In  state  A  -  a 
state  of  uniform  compression  -  the  membrane  shear  stresses  are 
aero;  in  state  B  they  are  finite  though  small,  Hence,  in  com¬ 
puting  the  membrane  shear  stress  contribution  to  JIq,  cognisance 


0 


a. 


should  be  taken  of  this  growth.  The  strain  energy  contribution 
will  be 


U 


s 


s: 


(shear  stress ^(engineering  shear  strain)^  dV. 


Again,  as  with  the  variation  in  the  normal  membrane  stresses 
in  going  from  state  A  to  state  B,  frequently  this  quantity  is  re¬ 
garded  as  small  in  comparison  with  Hg  and  is  neglected.  Never¬ 
theless,  we  will  retain  this  term. 

(c)  That  arising  from  bending  moments  and  torques  induced  in 
the  change  of  state. 

(d)  That  arising  from  the  work  done  by  the  external  force  system 
while  the  structure  deforms  from  the  "plane"  state  to  the  adjacent 
"bent"  state.  During  this  deformation  the  external  hydrostatic 
pressure  remains  constant  in  magnitude  (and  even  if  it  did  not 

it  would  be  assumed  that  it  did),  so  that  the  work- done  is  equal 
to  p  times  the  change  in  volume  occupied  by  the  shell. 

(e)  Strain  energy  in  the  reinforcing  rings. 


1*  Computation  of  potential  Energy  Contributions. 

Since  we  only  are  concerned  with  state  B  it  is  convenient 
to  take  stressed  state  A  as  datum  for  purposes  of  measurement  of 
the  additional  displacements  involved  in  going  to  the  "adjacent" 
state.  Accordingly,  we  assume  that  in  state  A  the  cylinder  is  of 
r’&dius  ft*  wall  thickness  h,  and  is  of  length  L  center  to  center 
of  ring  frames. 

bet  u, v,w  be  the  displacement  components  in  state  B  measured 
from  zero  in  state  A.  These  displacements  are  taken  as  positive 
according  to  the  scheme 


u  in  the  direction  of  x  increasing, 
v  in  the  direction  of  s  increasing, 
w  in  the  direction  of  r  increasing, 


where  the  coordinate  system  x,s,r  is  as  shown  in  Fig*  1.  As  can 
be  seen,  w  is  thus  positive  when  measured  in  the  direction  of 
the  inward  normal  to  the  shell  curved  surface. 

Let  the  energy  contributions  associated  with  (a),  (b),  (c), 
(d),  (e)  of  §2  be  Ua,  Ub,  Uc,  Ud,  Ug,  respectively.  Then,  with 
barred  quantities  denoting  the  stress  and  moment  resultants  (per 
unit  length  of  arc  on  the  middle  surface  of  shell  wall)  existing 
in  state  A,  and  unbarred  quantities  denoting  the  changes  in 
those  resultants  in  going  from  state  A  to  state  B,  we  have: 


(a):  Ua  «."mxxe]OC*Baaeaa)d4 


e  +N  e  )dA 
xxxx  ssss' 


(3) 


where 

flT  a  constant  , 
xx  ’ 

1  ft  *■  constant  , 
1  ss 


1 

J 


<u 


f"xx 


l 


Nxx^exx^  » 
Nss(ess)  > 


(5) 


and  the  region  of  integration  A  is  the  surface  area  of  the  cir¬ 
cumferential  ''middle  surface"  of  the  shell  wall. 


(b): 


Ub  ’  kS  . 


(6) 


where 


(12) 


p 


11. 

In  addition,  we  have  the  stress-strain  relations 
°*x  *  CE/tl-v2)J<exx  +  vess)  , 

“ss  "  CE/(l-va )3(esa  +  ve^)  , 

®xs  -  [E/2(l+v )]exs  . 

so  that  the  stress-resultants  are  given  by 

^xx  “  [E/(l-v2 ) ] (exx  +  vesg)  > 

Nas  -  [E/(l-v2)](ess  +  v&xx)  , 

Nxs  -  [Eh/2(l+v)]exs  . 

Substituting  (12)  into  (13)  then  gives  the  stress- resultants  in 
terms  of  displacements. 

In  the  second  integral  in  expression  (3)  we  normally  disre¬ 
gard  terms  of  degree  higher  than  the  second  so  as  to  obtain  a 
linear  problem  when  finding  approximate  solutions  by  means  of 
the  Rayleigh-Ritz  technique.  On  so  doing,  we  obtain,  therefore, 

®a  -  ■»  (Hp/2)  S  {ux  +  2vg  -  2(w/R)  +  (l/2)wx2  +  wg2}dA 
A 

+[Eh/2{l-v2)]  /  {u2+v|+(w/R)2-2vg(w/R)+2v(uxvg-ux(w/R)  )Jd£ 

A  (U) 

where  now  for  convenience  we  use  the  notation 


p 


12. 

For  we  have: 

“b  -  (1/2)  •£  *»•»<“  - 

<=  (l/2 ) [Eh/2 (1+v)]  s  (vig  +  vx  +  wxwg)2dA  ,  (15) 

A 

=  (1/2 ) [Eh/2 (l+v)]  /  (u.  +  v  ) 2dA  (15a) 

A  H  * 

on  neglecting  terms  of  degree  higher  than  two. 

To  oonpute  Uc  we  need  expressions  for  the  quantities  •  •  • , 
in  terms  of  displacement  components  u,v,w,  and  quantities  M^, 
in  terms  of  curvatures  or  displacements  via  stress-strain  rela¬ 
tions.  Since  these  terms  have  not  been  mentioned  in  the  discus¬ 
sion  by  Kennard  we  assume  the  results  of  PIBAL  number  167  without 

comment  or  elaboration  (PIBAL  number  I67,  page  4),  and  so  have 

Uc  =  (D/2)  f  {iXx  +  X/  -  2(l-v)(xxX  -  f2 )}dA  , 

where 

*X  8  wxx  ' 

^  ■  “ss  +  “/R2  • 
f  B  wxs  +  (vx/2R)  -  (us/2R)  . 

Hence 

Uc  -  [Eh3/?.4(l-v2)]/  {(w^  +  wsg  +  w/R2)2 

-  2(l-v  )[wxxwsg  ^  ww^/R2  -  (wxs  +  vx/2R  -  ua/2R)ajjdA. 

(16) 


To  compute  the  contributions  from  the  end  forces  and 
from  the  radial  forces  will  he  considered  separately.  We  note 
that  the  end  thrust  produces  (among  other  things)  a  stress  re¬ 
sultant  Hxx  of  magnitude  -  Kp/2.  If,  now,  we  consider  an  ele¬ 
ment  of  the  curved  surface  of  the  shell,  of  length  dx  add  (circum¬ 
ferential)  width  ds,  and  defined  longitudinally  by  the  two  points 
AB  as  shown  in  Fig.  2,  we  have: 


displacement  in  x  direction  of  point  A  ■  u^  , 
displacement  in  x  direction  of  point  B  *  u^+du^  . 

Hence,  for  the  element  in  question,  the  work  done  by  If  is  given 
by 

AUde  3  (RP/2)tuA  **  (uA+duA^ds  > 

and  so 

Ude  "  -  (RP/2)  S  VA  ‘  (17) 

A 

In  passing,  and  although  not  strictly  relevent  here,  we  note 
that  if  the  cylinder  were  short  and  closed  at  the  ends  the  work 
done  by  the  end  thrust  would  be 


P  S  [u]  ■  dA  . 
end  X“ 


For  the  second  contribution  to  Uj  we  note  that  if  any 
normal  force  p  at  a  given  point  on  the  curved  surface  of  the 
shell  moved  with  the  point  during  the  displacement  of  the  point 
on  taking  up  state  B,  then  (and  as  is  usually  done  in  Elasticity) 
the  work  done  by  such  a  force  at  that  point  would  be,  simply, 


p 


I 


p 
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Following  Kenuard*,  however,  we  consider  this  matter  further. 


£.  Work  done  by  Lateral  Hydrostatic  Pressure. 


In  more  detailj  we  consider  a  point  0  located  on  a  lateral 

surface  generator  in  state  A,  and  displaced  to  a  position  D  in 

state  B,  Fig.  3«  Bet  D*  be  the  projection  of  D  on  the  generator 

passing  through  G.  .Then,  if  the  normal  applied  force  moved  with 

C,  its  normal  displacement  would  be  w  .  If,  however,  the  force 

c 

originally  at  C  remained  in  that  same  position  and  did  not  move 
with  G,  its  normal  displacement  would  be  [w]  .  This  we  can  take 

V 

to  be  the  situation  as  far  as  the  lateral  surface  hydrostatic 
pressure  is  concerned.  From  Fig.  3,  on  allowing  for  the  circum¬ 
ferential  displacement  vc  of  point  C  (not  shown  on  the  figure), 
we  see  that  to  within  a  first  derivative  correction  [w]c  is  given 
by 


j>]c  -  wc  -  uc  (wx)c  vc  (ws)c  . 


Using  this  result  we  can  compute  the  change  in  volume  of  a 

cylinder  of  original  radius  R  and  length  I>.  The  original  volume 

is  IIR2dx  ,  The  final  volume  is 
o 

L  2D  p  , 

V  =  S  |(l/2)  S  [R  -  (w  -  uwx  -  vws)]  d<pjdx  , 


or 


I,  .  2  HR 

V  =  S  |{l/2R)  /  [R  -  (w-  uwx  -  vws)]2dsldx  . 

O  0 


*  Kennard’s  discussion  on  fIBAL  number  167  is  included  for  con¬ 
venience  in  this  reoort  as  Appendix  IV. 
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P 


15. 

Thus  the  change  in  volume  is  given  by 

V  -  /{w  -  (w2/^)  -  uwx  -  vwJ-dA 
A 

where  terms  of  degree  higher  than  the  second  have  been  neglected. 
Hence  we  see  that  the  work  done  by  the  radial  forces  is 

Udr  “  P/(w  "  (w2/2R)  -  uwx  -  vwg|dA  ,  (Id) 

A 

where,  as  usual,  A  is  the  area  of  lateral  or  curved  (middle)  sur¬ 
face  of  the  shell  in  state  A. 

To  all  the  foregoing  we  should,  for  completeness,  add  energy 
terms  due  to  strain  in  the  end  reinforcing  rings  frames.  This 
we  will  not  do  here,  however;  and,  when  such  terms  are  required 
we  will  take  them  direct  from  PIBAL  number  16?. 

YJith  the  above  material  as  framework  we  are  now  in  a  position 
to  discuss  material  contained  in  PIBAL  number  167  together  with 
the  remarks  on  that  material  by  Dr.  Kennard. 
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PART  2.  DISCUSSION 

6.  Quantity  (Ua-U^)  derived  here  and  by  Kennard. 

In  his  work  Dr.  Kennard  obtained  some  energy  quantities 
which  differ  from  those  given  in  PIBAL  number  167.  Essentially, 
the'  differences  occur  in  the  energy  contributions  Ua  and  U^.  We 
consider,  then,  the  quantity  Ua-Ud. 

From  (14),  (17),  and  (16)  we  have 


Ua-Ud  »  -  ( Rp/2 )  «/’{ux+  2vg-  { 2w/R ) + ( wx2/ 2 )+ wg a} d A 
~~  A 

+£Eh/2 ( 1-v2 ) ]  /{u2+v|+ ( wz/R2 ) - ( 2 vaw/R )+2v (uxvs-uxw/R ) }dA 
Pi 

+(Rp/2)  J*  uxdA-(Hp/2)  /{(2 w/R)^(w2/R2)^(2uwx/R)-(2vws/R)}dA. 
A  A 


Further,  we  note  that  v  is  periodic  in  the  circumferential  direc¬ 
tion,  with  a  consequence  that  J*  vsdA  has  magnitude  aero.  Thus 

A 

we  have 


ua“ud  ”  RP  u\r(w2/2R£1)+(uwx/R)+(vws/R)-(w2/4)-(^/2)JdA+(AUj  )b 
A 

where 

(AU?  »  [Eh/2(l-\^]  /[u|+v|+(w2/R2)-(2vsw/R)+2vUxvs-(2vuxw/R)]dA{ 

(19) 

This  is  the  result  obtained  from  the  foregoing  considerations 

only. 

On  making  obvious  changes  in  notation,  expression  (19)  is 
in  complete  agreement  with  the  corresponding  expression  (e)  given 
by  Kennard, 


p 
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7.  quantity  (Ua-Ud)  obtained  from  PIBAL  number  167 . 

After  a  little  manipulation  the  quantity  (Ua^d^  obtained 
from  the  work  in  PIEAL  number  167  is  given  by 

Ua-Ud  -  Rp  ur{.-<w/R)+(w2/2R2)+(w  wsg/2)-(w  ux/2R) 

A 

-  (wx2/4)-(vx2/4)-(ux2/4)}dA  +  ( AU  )fe  .  (20) 

Note  that  as  written  here  this  is  not  the  expression  (e*) 
quoted  by  Kennard  in  his  work  as  the  PIBAL  result  in  question. 

It  should  be  remarked  here  that  in  (19)  no  account  has  been 
taken  of  the  membrane  3hear  stress  contribution  U.  .  On  the  other 

D 

hand,  as  obtained  from  PIBAL  number  16?  as  will  be  shown  shortly, 
(20)  comes  from  (but  does  not  include)  quantities  which  take  into 
consideration  such  effects. 

From  a  comparison  of  (19)  and  (20)  there  are  obvious  dif¬ 
ferences,  and  it  would  appear  that  the  latter  expression  is  in¬ 
correct.  Nevertheless,  although  a  portion  of  the  theoretical 
treatment  presented  in  PIBAL  number  167  seems  in  error,  it  will 
be  shown  later  that  the  results  and  conclusions  given  in  that 
report  remain  essentially  the  same  as  those  obtained  from  (19). 
Before  doing  this  it  is  of  interest  to  discuss  the  deviation  of 
(20)  further. 

The  PIBAL  number  167  equivalent  of  (20)  is  made  up  from 
equations*  (3),  (25)  and  (25)  of  that  report.  For  completeness 


*  The  numbers  of  the  PIBAL  report  equations  are  here  shown  as 
berred  numbers  for  clarity.  In  that  report  the  numbers  are  the 
same  but  they  are  not  barred. 
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they  are  repeated  here.  We  have 

2n 

(i)  U  -  [Eh/2(1- v2)]  i2  +U  -w)2  +  2  y  u  (v  ~w)  +  (1-  v)6r+u  ff  2]  d£  dci 

6  o  0  5  <P  £  <P  l  q>  J 

(3) 

where 

U  is  "the  energy  of  stretching", 

^  m  x/a  or,  in  our  notation,  -  x/R, 

Lq  -  L/R; 

271  L_ 

(ii)  V  -  ~  (Ph/2)/  S  (u«  +  v*  +  w*  )  d£  d9  ,  (20) 

l  ^  Q  S  £  5 

where 

V  is  the  "potential  caused  by  axial  pressure",  and,  on  interpretation 
of  quantities  given  on  Pibal  pages  17  and  19, 

P  *  pR/2h, 

and  on  page  9  of  the  Pibal  report  P  is  defined  as  the  "axial  pressure"! 

2n  L  4  _ 

(iii)  1?2  -  U/2^*  S  pw  {2  +  t-  (w  +  w^  )/R>-  R  rtx  09  ,  (25) 

where 

H  is  the  "v.ork  done  by  the  radial  pressure." 

We  note  that  (3)  as  ,iust  presented  and  as  used  in  Pibal  number  167 
represents  the  total  strain  energy  due  to  all  membrane  stresses,  i.e.  shear 
es  well  as  normal,  end  we  see  that  7^  as  defined  abov<>  is  the  negative  of 
the  work  done  by  the  pressure  acting  on  the  ends  of  the  cylinder. 


p 
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In  virture  of  the  definition  of  Uq  .just  given,  it  follows  that  UQ 
reoressrrtc  out*  expression  (U  +  U,  ).  In  consequence  of  this,  and  the  fact 
that  Ufe  is  not  present  in  (19),  it  is  convenient  to  write  Ue  as 

”e“Ue  +Ue  > 

1  2 

where 

UA  -  (Bh/2(jU>  v2  )3  S™  f°  f uZ  +  (v  -  w)2  +  2vu  (v  d?  , 

“  o  o  u  5  9  t  J 

2ft  I*  2  1 

U  -  CSV2(1-V2))/  £  f  (1**  V  )(v  ♦  U  )  /2}  d?  dep  . 
e  oo'  F  9 

2  ^ 

To  transform  these  to  our  notation  we  use  the  substitutions 


from  which 


K  -  x/R,  a  «  Rq>  , 


d£  =«  dx/R,  dq?  »  ds/R  , 


and,  on  noting  that  Lq  -  L/R,  it  follows  that 


U  -  [£h/4(l+  v  )3  Z}R  S  (Rv  +  R»  J2  ds/R2  , 

e-  o  o  *  a 


’nK  L 

[Eh/4(1+  v)3  £  S  (v  +  O2  dr.  de  , 

0  0  X  B 


Uv 


Consequently,  in  obtaining  (?0) ,  i.o.  the  equivalent  of  (Ua  -  U^), 


we  must  udh  U  instead  of  II  , 
*1  6 


3 


P 
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Using  the  same  transforms lions  for  l  end  cn  m  obtain  for  U  the 


expression 


J  JjTfR  i*  „  -  _  Z 

u  »  [i'.h/2(i-.v  )]  s  [fr  u*  +  (Rv  -  ft)"  +  2v  Ru  (Rv  -  vt)  Vdx  ds/R  , 

«  o  o  *  3  xej 

2  ZnR  L 

a  [£h/2(l^v  )]  J"  v,'  fu£  +  \r -( 2wvs/R)+ ( wa/R?' )  +  2  vuxvs-?.v  wu/R  j-  dx  do  . 

(21) 

In  similar  manner  we  can  compute  ,  noting  that  we  require  this 

difference  since  V  is  the  negative  of  the  work  done  by  the  "end  force" » 

'Ha  have 

v  -  w  3  -  (Ph/2)}71  f°  (»*  +  v*  +  w*  )  dr  dtp 
*  00  t,  n 

*  -  U/2)f  /  P»i  {2  +  a  *-  (w  +  w  )/R  }  R  dx  dtp  , 

rt  n  v  X  J)  W  J 


o  o 


and,  on  noting  that 


this  becomes 


Ph/2  -  (PR/2h)(h/2)  «-p  R/4  , 


2nR  L 


V  ~  «  -  (pH/4)/  /  {(R  u/  +  (R  v/  +  (R  wx)Z  J  dx  ds/R‘ 

2nR  L 

«-  (l/2 )  J‘  /  pw  {2  +  -  {w  +  RZvif.a)/R }  fix  ds, 

awR  L 


"  -  PR//  (  [(^4MvJ/*Mi^4)}  +  *{2  mix~-(^^»ss)/R} /2R& 

*3®  L  ^  22  .  ^ 

P  Rj;  ry  {  ~(\/ ^)“(v/-0^x/4)-(w/SH^/21t)+  (*• /2R  j*  (r.wsg/2)}d  ' 


(22) 


p 
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Thus,  on  comparison  we  see  finally  that  (21)  is  the  bame  as 

( A’J  )b  of  (20)  and  (22;  i&  the  3ama  3s  the  Quantity  under  the  integral 

sign  in  (20),  so  that  the  quantity  from  Pibal  number  167  to  he  compared 

with  (tia~Ud) ,  n-arraly  L*e  *■  V  -  W  ,  is  given  by  (2C)  as  was  stated, 

112 

3.  /Motivation  of  durivation  in  Pibal  number  167. 

As  the  authors  of  Pibal  report  number  167  are  now  not  available,  it 
cannot  he  stated  with  any  certainty  what  were  the  bases  on  which  tho 
expressions  contributing  to  (20)  were  darived. 

It  is  conjectured  that  for  U  of  Pibal  number  167,  that  is  the  energy 
contributions  from  all  membrane  stresses,  only  the  second  integral  term 
in  (3)  was  used  instead  of  the  full  two-term  expression  of  (3);  in  addition, 
of  course,  the  quantity  Ub  as  given  by  (6)  was  also  used.  The  second  term 
of  (3)  rerpesenls  energy  due  only  to  the  change  in  magnitude  of  tho  normal 
membrane  stresses  in  going  from  state  A  to  state  D,  and  the  If  of  Pibal 
number  167  apparently  neglects  the  existence  of  normal  membrane  stresses 
already  present  In  state  A, 

In  the  expressions  fox-  the  work  done  by  the  end  forces,  the  Pibal 
number  167  expressions  can  bo  derived  by  assuming  that  t,ha  complete  non¬ 
linear  strain  expression,  nuusoly 


j  •  U  f 

XX  (X. 


(•4* 


X 


w2)/2 

X ' 


holds,  and  then  by  imposing  tho  condition  that  the  buckling  process  is 
essentially  inextsnsional,  i,e.  a  0,  etc.  It  is  usually  considered 


p 
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that  the  noiv»Itaear  terms  u^,  vx  sru  of  higher  order  than  »„  ,  and  in 

the  Pibal  report  all  tern's  arising  from  these  were  subsequently  eliminated 
from  the  final  results.  That  this  was  the  manner  of  derivation  in  the  Pibal 
work  is  again  pure  conjecture, 

The  Pibal  expression  for  the  work  done  by  the  radial  forces,  i,o.  (25), 
is  open  to  criticism  in  that  it  seems  that  to  compute  the  local  change  in 
radius  of  the  original  circular  cylinder  the  radius  of  curvature  of  the 
deformed  cylinder  (state  B)  was  used.  If  this  was  so,  it  is  incorrect. 

Effect  of  errors  in  original  derivation. 

Having  obtained  (20),  and  discussed  the  possible  manner  of  the  original 
derivation  of  that  expression,  we  now  consider  the  difference,  D  say,  be¬ 
tween  (20)  and  (19).  '/riting 

B  »  (19)  -  (20) 

we  have 

D  »  Bp«"C  {(w  \xJZR)  +  (v^AMu^/4)  -  (w  wgg/2)  +  (w/B)} 

+  {(u  wx/ft)  +  (v  wg/R)  -  (w^/2)  }  ]<JA  , 

whore  the  quantity  in  the  first  set  of  curly  braces  consists  of  the  terms 
left  from  Pibal  number  16",  and  that  in  the  second  set  of  curly  braces  is 
from  Kannard's  work  when  terms  common  to  (19)  and  (20)  are  cancelled. 

As  will  be  discussed  in  detail  later,  in  the  first  quantity  certain 
individual  terms  are  subsequently  neglected  in  Pibal  number  167  on  the 


ba6 


In 

in 

The 


wen 


whej 


and 
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basis  of  a  smallness  argument.  The  terms  so  neglected  are 


r-  ux/2R  ,  v^/4  ,  u^/4 

In  a  similar  manner,  that  is  for  comparable  accuracy  (again  to  be  discussed 
in  detail  later)  two  terns  in  tho  second  quantity  can  also  be  neglected. 


These  are 


u  i»x/R  ,  v  wg/R  . 


Thus,  essentially,  the  difference  between  (19)  and  (20)  is  given  by 


D  -  -  (Rp/2)  f  *  -  2w/R  6A  .  (23) 

A  as  s 

From  Pibal  number  167  the  chosen  displacement  coefficient  functions 


u  *  A  cos  m<J>  cos  -  A  cos  (ms/R)  cos  (  \  x/R)  , 


v  «  B  sin  mcp  sin  If.  a  B  sin  (ms/R)  sin  (  A  x/R)  , 

w  «  C  cos  mcp  sinAg  »  C  cos  (ms/R)  sin  (Ax/R) 


where  m  is  the  number  of  curoumferential  halfwaves. 


A  m  nnB/L  , 


and  n  is  the  nuubor  of  half-waves  i.n  the  axial  direct-xe. 

Substituting  (27)  in  (23)  w a  obtain  thu  result,  mentioned  earlier. 


D«-  (Rp/2)  S  ft*  wss  +  vT  -  2w/r}  dA  «  0 
A 


(23a) 


This  is  discussed  in  more  detail  in  Appendix  ITX. 
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Thus  we  conclude  that,  although  in  no  way  a  justification  for  the 
theoretical  treatment  given  in  Pibal  number  16?,  and  depending  on  a  smallness 
argument  still  to  be  presented,  the  simplified  formula  (43)  given  in  the 
Pibal  report  stands.  The  formula  is 

X^l-v^/Cm2 +  X.^  +  Rc  +  k(m2 +\2) 

$  , - - - >  (43) 

(i!?-l  +  X2/2) 

where  tho  various  quantities  and  symbols  concernedwill  be  mentioned. in 
detail  later. 

It  can  further  bo  stated  that  the  numerical  work  relating  to  Pibal 
number  167  was  based  on  eouatp-on  (43) . 

The  energy  terras  of  Pibal  number  167  were  used  in  other  Pibal  reports, 
namely  Pibal  numbers  I69,  177,  182  and  189.  These  reports  extend  tho 
previous  work  by  considering  other  shell  constraints  at  the  frames,  by 
discussing  various  aspects  of  shell  behavior  and  (in  189)  by  considering 

t 

the  overall  buckling  of  a  reinforced  shell  oetween  bulkheads.  In  all 
reports  deflection  functions  are  used  which  satisfy  equation  (23a). 

Smallness  arguments  similar  to  those  presented  here  are  used  in  these 
reports  in  the  derivation  of  the  simplified  formulas  7;hich  were  used  for 
numerical  calculations.  On  this  basis  it  seems  that  the  working  formulas 
and  numerical  work  of  the  Pibal  reports  using  the  energy  expressions  of 
Pibal  number  167  are  valid. 


p 
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PART  3.  DETECTIVE  WORK 

10.  Essential  stops  in  obtain:  :.e  (43)  of  Plbal  number  167. 

Having  shown  that,  on  neglecting  certain  terms,  the  main  result 
obtained  in  Pibal  number  167,  namely  (48),  still  holds,  we  now  discuss 
the  justification  of  the  approximations  resorted  tp.  To  do  so  it  is 
necessary  to  outline  the  work  presented  in  the  Pibal  report. 

In  the  main  body  of  the  present  report  it  has  not  been  emphasized  that 
the  cylinder  being  considered  is  reinforced  with  ring  frames  and,  in  fact, 
expressions  for  the  strain  energy  in  the  rings  have  not  yet  been  given. 

These  we  will  not  discuss  here.  We  present,  instead,  the  complete  expres¬ 
sion  for  the  total  potential  energy  as  used  in  Pibal  number  167.  Note  that 
this  includes  the  supposedly  incorrect  contribution  (20).  Consequently, 
after  discussing  the  methodology  and  manipulation  used  in  the  Pibal  report, 
we  will  than  draw  corresponding  conclusions  obtained  on  using  (19)  instead  of 
(20). 

From  Pibal  number  167  wo  have  for  the  total  potential  onergy  ti  a  expres¬ 
sion 

°T  -  "c  ‘  Ub  *  W  V,  (55) 

where  the  individual  contributions  Uy,....,  are  given  by  (3),  (4),  (13), 
(20)  and  (26)  respectively.  Here 

UT  »  total  "potential",  i.e.  total  potential  onergy, 

Ue  «  "strain  energy  of  stretching",  i.e.  membrane  stresses  strain  energy, 

*  "strain  energy  of  bending", 

Uy  -  "energy  stored  in  a  ring",, 

V1  •  "potential  caused  by  the  axial  pressure", 

V  m  "potential  energy  due  to  the  radial  pressure". 
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In  Pibal  number  167  the  work  then  follows  tho  usual  motivation  of 
the  Rayleigh-Ritz  process  for  computing  approximations  for  eigenvalues. 

The  chosen  displacement  functions  (27)  are  inserted  in  (34)  and  the  indi¬ 
cated  integration  carried  out. 

Conditions  necessary  for  rendering  U^,  stationary,  namely 

sut/3a«  auT/aB-  au^c-o  ,  (35) 

are  then  imposed.  Those  lead  to  a  set  of  three  simultaneous,  homogeneous, 
algebraic  equations  linear  in  A,  B,  C.  For  a  non-trivial  solution  to  exist 
the  familiar  condition  that  the  determinant  of  tho  coefficients  of  A,  B,  C 
must  bo  zero  is  than  used,  which  loads  to  the  eigenvalue  detorminantal  equa¬ 
tion  (39),  namely 

fx2  +  m?(l~v)(l*k)(l/2)  {- \u£ (l/2)(l+vWk/2)(l-v))}  (Xv-(l-v)  Xkm2 

*  "A  "  X*®2!  *SAC**,X/2) 

{  -  X  m[  ( 1/2)  (  1+v )-( k/2) (1-v) ] j  fn^l-v)(l+k)(X2/2)-£«>  \  f-m(l*(l-v)X2  k^j 

2* 

(Xv  -(1-v)  X  km2  {-fl£  1+  ( 1— v)  X“k]'j  fl+kC  (m2+XZ)  +l-2m2  -2vXZ] 

+  V  *1  k/Z^  +  Rc 

(39) 

It  is  then  shown  in  Pibal  number  167  that,  for  a  cylinder  for  which 
the  distance  between  adjacent  ring  frames  becomes  largo  (-►<»),  the  solution 


p 
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of  (39)  yields  the  classic  Bryan  solution  ior  an  unstiff enod  cylinder* 

p  =  (ni2  -1)  £iih3/l2R3(l-.v2)|-  .  (44) 

In  this  treatment  various  quantities  were  introduced,  in  particular 
we  have 

K  -  Eh/(l-v2)  ,  k  «  hz/l2R2  |  (29) 

«  .  p  R/K  (  -  pR(l  -v2)/Sn)  , 

$  -  Ph A  (  -  pR(l-V2)/ZEh) 

for  the  hydrostatic  loading  case  discussed  bar*,  so  that 

«  /2  .  .  (43) 


As  mentioned  in  Pibal  number  167,  "the  complete  expansion  of  the 

determinant  (39)  leads  to  a  cubic  oquation  in  $  and  $  Consequently 

1  2 

as  a  first  approximation  (see  Appendix  1),  in  the  expansion  powers  and 
products  of  and  k  greater  than  tbo  first  can  be  neglected. 

This  leads  to  tho  approximate  expansion 


or,  using  (43), 


C  +  C  k  ■  C  $  +  C.  $ 

1  2  3  1  *  2 


Cj  +  Cg  k  ■  (C3  +  C^/2) 


(42) 


(44) 


Expressions  for  ^  3  are  given  by  (£4a),...,  (44d)  of  Pibal 

number  167.  On  inspection  of  these  it  is  apparent  that  for  m  large  (e.g.  m>  8) 
each  expression  contains  one  terra  dominant  in  comparison  with  the  remaining 
terms  of  that  expression.  This  dominant  term  is  independent  of  the  ring 


parameters.  In  the  expression  for  there  also  occurs  a  second 

dominant  term  of  similar  magnitude  but  involving  the  ring  parameter  R  • 

c 

If  only  the  dominant  terras  are  retained  and  used  in  (44)  obtain, 
finally,  the  approximate  result  (48).  Further,  if  in  this  equation  the 
terra  containing  the  ring  parameter  is  deleted,  the  simplified  von  Hises 
result  follows,  (45)*  Actually  there  is  a  small  discrepancy  here,  for 
in  lieu  of  the  denominator  tom  (ra2  -1  +  (X2/ 2)}  of  (45),  the  von  Mice1  s 
result  neglects  the  -1  contribution  in  comparison  with  m  . 

11,  Justification  or  neglect  of  certain  small  Pibal  toms. 

The  foregoing  discussion  retraces  the  main  steps  taken  in  Pibal  number 
167  in  obtaining  (48).  In  the  process  many  small  quantities  were  neglected. 
Among  them  are  terms  related  to  the  terms  mentioned  earlier,  i.e. 

*  ux/2R  ,  v2  /4  ,  u2/4 

We  now  discuss  these  three  terms  in  more  detail. 

If  (39)  is  written  as 

iaiji  "  0  »  i»J  *  1iZ-3  » 

it  is  easy  to  trace  the  following  connections* 

(i)  In  a.,  the  quantity  \2f.  comes  directly  from  the  term  u2  of 

11  2  £ 

(20),  that  is  from  . 

(ii)  In  a„_  the  quantity  \2$  comes  directly  from  the  torn  v2  of 

22  2  £ 

(20),  that  is  from  v*  . 

(Hi)  In  a  1  3  and  a  the  quantity  /Z  coitus  directly  from  the  term 


p 


p 
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Further,  it  is  simple  to  show  that  the  H  term  drops  out  from  the 

denominators  of  the  terms  u2  .  v2  • 

x  x 

In  the  expansion  of  (39)  these  three  quantities  lead  to  terms  in 

$  and  $  of  degree  higher  than  the  first  (and  so  can  bo  neglected) 

and  to  contributions  in  the  non-dominant  tertns  of  C  and  C  .  They 

3  4 

make  no  contribution  to  the  dominant  terms  retained  in  deriving  (48). 
That  this  is  so  is  shown  in  detail  in  Appendix  I. 

12j  Eigenvalue  Determinant  for  the  Ksnnard  Energy  terms. 


Having  discussed  the  technique  and  reasoning  used  in  Pibal  nunfcur 
167,  we  non  apply  the  same  thoughts  and  degree  of  approximation  to  the 
onorgy  quantities  deduced  by  Kennard. 

For  the  total  potential  energy  wo  have 


where 


Ua  + 


+  -  U .  + 

bed 


U 


» 


is  the  energy  due  to  the  normal  membrane  stresses,  (14), 

is  the  onorgy  due  to  the  shear  membrane  stresses,  (15a), 

U  is  the  energy  due  to  the  moments,  (16), 
c 

Ud  is  the  work  done  by  the  external  force  system,  (17)  +  (18), 
Uq  is  the  energy  due  to  ring  deformation,  (13) . • 


Thus,  in  this  0^  we  use  expressions  from  Pibal  number  167  for  Uc  and 
Uq  ,  after  transforming  ttoio  to  x,s  coordinates. 


p 
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p 


Oil  following  the  procedure  proscribed  in §10,  in  lieu  of  (39)  we  obtain 
tho  determinants!  equation 

fx2+  ma(l~v)(l+k)/2  -Xa  {(l+v)-k(l-v)  }  /2  (Xv  -kXm2(:h-V) 

♦  Ba]  *  *1C  *  *•  •,)  .  . 

-  Xm  {(l+v)-k(l**-  v )  J  /2  fm2  +  (l-V)(l»k)X2/2]  hj£1+  ®,  +k  X2(!Uv)] 

f-  Xv-k  Xm^l-v)  -nC  1*#  +kX2(l-V)3  (l+k  ;  {(m2  +X2)  2  +l_2m2  -  2  v\ 

j 

♦  Rac  ♦X®1  ]  +  Rc  -  (m2  **3*  X*/2)®,} 

(24) 


The  connate  expansion  of  (24)  leads  to  a  quadratic  equation  in  $ 
and  $  ,  where  those  quantities  are  tho  same  as  previously  defined* 

M 

Again,  as  a  first  approximation,  powers  and  products  of  $  and  k  can  be 
neglected  in  the  expansion  and,  so  doing,  we  obtain 

f  (25) 

Here  C  and  Cfi  are  new  coefficients  corresponding  to  C1  and  C  2  of  Pibal 

nuabor  167,  and  the  now  coefficient  C? takes  tho  place  of  the  previous 

(C3+  C4/2).  Expressions  for  C5  ,  C6  ,  Cy  are  given  in  Appendix  2* 

As  before,  on  inspection  of  these  expressions  it  is  apparent  that  for 

m  large  (>8,  say)  each  one  contains  one  dominant  term  which  is  independent 

of  the  ring  para  outers,  and  in  there  ia  also  a  second  term  of  similar 

magnitude  involving  the  ring  paramo tor  R  .  Further,  it  can  be  seen  that  ‘ 

c 

those  dominant  terns  are  the  same  03  those  in  the  expressions  for  C  ,«.*,  C  • 

1  4 

Again,  on  retaining  only  tho  dominant  toms  in  C  ,  C  ,  C  ,  (25)  loads 

5  6  7 

directly  to  tho  previous  result, ( 48) * 


p 


p 
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13,  Justification  for  neglect  of  certain  siriall  Kannard  t6nas» 

In  the  immediately  preceding  derivati  on  of  (48)  cany  quantities 
have  boon  neglected.  Among  those  so  treated  are  quantities  arising 
from  the  terms 

u  nx/R  ,  v  w^R  , 

As  before,  we  nos  discuss  these  tso  terms  in  detail. 

If  so  srite  (24)  as 

lbijl  "°  *  -  1»2»3» 

it  is  a  simple  matter  to  trace  the  following  connections* 

(i)  In  b  and  b  the  quantity  combs  directly  from  tho  tom  u  s  /R  ; 

13  31  1  x 

(ii)  In  b  and  b  tho  quantity  m  i  coiosa  directly  from  the  term  v  w  /R  . 

23  32  1  » 

In  the  expansion  of  (24)  these  quantities  lead  to  second  degree  terms 

in  <6  and  4  (and  so  can  be  neglected),  and  to  contributions  to  th6  non- 
I  2 

dominant  toms  in  C  *  They  make  no  contribution  to  the  dominant  terms 
retained  in  deriving  (48).  This  is  shown  in  Appendix  II. 

14.  The  Kennard  Determinant,  and  the  Bryan  Solution. 

Vfe  consider,  now,  tho  effect  on  tho  Kennard  dwtonainanfcal  equation 
(24)  of  allowing  the  distance  between  the  ring  frames  to  become  large. 

On  lotting  L  -->00  or,  what  is  the  same  thing,  X  ->-0  in  (24)  we  obtain 

m2  £  1  ^  (1  +  ^  J2  +  $  (l~ra2)  +  k(m2  -1  )2|  *  0  , 

(26) 


p 
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p 


a 

If  we  expand  this,  neglect  4^  in  comparison  vdth  *  ,  end  substitute 
for  <15  in  terms  of  p,  we  obtain 

p  =<  [nt2  -  3  +  4/(m2  +1)}  (  Eh3/12R3(i~  v2)  j  , 

which  is  not  the  Bryan  solution  (41).  If,  however,  we  write  simply 

(1  +  fc^si 

in  (26),  we  do  obtain  the  Bryan  solution.  Obviously  in  so  doing  we  are 
not  approximating  consistently. 

The  torn  producing  this  inconsistency  can  easily  be  traced  back  to 
earlier  considerations.  In  (26)  the  term  (1  *  )2  conns  from  the  terms 

^23»  ^32^°  (24)  on  letting  \->»0.  And  in  b23  ,  b3JJ  the  quantity  $  in 
the  term  (l  +  *1^  )  comes  directly  from  the  term  v  w ^/H  which  was  intro¬ 
duced  originally  in  the  expression  for  the  work  done  by  the  lateral  thrust. 
The  other  term  of  similar  character  introduced  in  t'net  samo  expression, 
namely  u  w^/R  dcos  not  survive  here  since  its  end  product  in  (24),  namely 
,  van! shoe  with  \  . 

15.  Further  Discussion  of  7«ork  Done  by  internal  Forces. 

The  computation  of  the  work  done  by  external  forces  warrants  further 
discussion.  Tie  have  just  seen  that  use  of  the  expression 

Udr  "  p  f  w  **  w2/2R  “  u  \  ~  v  VJ'M  •  (IS)  bis 

A 

for  the  work  done  by  the  forces  normal  to  the  lateral  curved  surface  of 
tha  shell  leads,  apparently,  tc  an  inconsitoncy  in  comparison  with  a  result 


p 


p 
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obtained  by  Bryan.  The  .inconsistency  arises  not  from  an  error  but  rather 
from  an  order  of  magnitude  disagreement-. 

The  tonus  u  w  ,  v  w  refer  the  final  position,  after  deformation  to 
state  B,  of  a  point  on  the  lateral  surface  back  to  its  original  position, 
and  so  allow  for  a  more  correct  integration  than  is  usually  the  case 
(in  so  far  that  a  slope  correction  has  been  added)  with  respect  to  x  over 
the  range  0<x<  L  which  defined  the  cylinder  in  state  A.  Nevertheless, 
even  when  using  non-linear  strain  components,  this  correction  is  rarely 
used  in  Elasticity  problems.  Usually,  even  though  deformation  lias  taken 
place,  no  allowance  is  made  for  the  change  in  area  of  an  element  of  sure 
face  and  integration  is  carried  out  over  regions  defined  by  original 
unstrained  boundary  positions. 

The  Potential  Energy  Theorem,  used  in  Elasticity  arid  here,  states  that 

* 

for  an  equilibrium  state  thu  potential  energy  U 

U  "  EG  ~  ^  Ti  ui  -ds  »  1  **  h  2,  3,  (27) 

SB 

is  stationary  for  arbitrary  variation  in  displacements.  Here 
- 1  {  °«  *«  ’  i>J  ■  1'2’3 


^ith 


Ji.j 


(•.  .) 
3J 


via  the  stress-strain  law,  and  T  ,  u 

.1  i 


i  »  1,2,3,  are  the  components  of 


*  For  the  present  we  ref or  all  quantities  to  an  orthogonal  Cartesian 
coordinate  systemj  sec  Sokolnikoff  "Theory  of  Elasticity*1 . 


p 
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surface*  traction  and  of  tho  corresponding  surface  displacement.  The 
portion  of  the  surface  designated  as  Sg  is  that  portion  for  which 

p 

a  6Ti  dS  =»  0  . 


Xh  particular  if  surface  tractions  are  specified  over  the  entire  surface 

of  the  body  in  ouost.ionj  then  5T,  s  0,  and  S_  <s  entire  surface, 

1  B 

In  these  various  expressions  the  usual  summation  convention  is 
understood,  and  in  all  applications  the  limits  of  integration  ore  the 
original  dimensions  of  tho  unstressed  elastic  body. 

If  wo  formally  apply  theorem  (27)  to  the  present  problem  the  term 
corresponding  to  the  "work  done  by  the  hydrostatic  pressure"  on  the  lateral 
surface  of  the  shell  is 

U'  -  /  pwdA  (28) 

dr  A 

whore,  as  usual  throughout  this  papor,  dA  *  ds  dx,  and  A  refers  to  tho 

curved  surface  area  of  tho  shell  in  state  A. 

An  expression  such  as  the  last  represents  a  computation  correct  only 

to  the  first  order  and,  truly,  the  Kennard  correction  turns  added  to  the 

normal  displacement  w  yields  a  quantity  U .  which  is  more,  accurate  than 

dr 

is  (28). 

As  an  alternative  to  the  Kennard  computation  one  could,  of  course, 
compute  the  work  done  by  the  lu„eral  throust  by  using  the  formula 


ti 

U. 

ar 


dA 


where  now,  on  writing  dA  =*  dx  ds,  the  range-  of  definition  of  x  would  be 
0  s  X  <  L  —  Al>,  where  AL  is  the  shortening  of  the  cylinder  on  going  from 
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state  A  to  state  V,  and  is  given  by 

AI»  «  u  -  u 
L  o 

However,  if  the  quantity-  U  *!  were  used  in  this  connection,  then  similar 

dr 

limits  of  integration  should  also  bo  used  in  all  other  computations  throughout 
this  work.  This,  of  course,  has  not  been  done. 

It  seems,  thon,  that  on  retaining  the  terras 


equi 

the 

cyl 

to 

-  equi 
von 
can; 


in  tho  expression  for  U^,  computation  throughout  the  treatment  is  not  in  ■ 

consistent  in  accuracy.  fio 

Again,  in  retaining  tho  limitB  0  <  a  <  2nR,  the  term  ii2/2^  in  (18)  is  ^ 

refora  to  the  change  in  the  length  ds  of  the  element  on  going  from  state  A  solj 

to  state  B.  Thus,  for  accuracy  consistent  with  that  in  all  other  surface 

area  computations  this  term,  likewise,  should  not  bo  used,  obti 

16.  Consequence  of  ds  correction  term.  witl 

,  2,  d0n' 

Since  in  all  the  foregoing  work  the  term  w  /2R  has  been  retained  it 

is  of  interest  to  trace  its  consequences,  and  so  discover  the  effect  of 

ignoring  it. 


On  substituting  the  deflection  expressions  (27)  in  the  total  potential 
energy  expression  (34)  and  performing  the  differentiation  with  respect  to 
the  arbitrary  parameters  A,  B  and  C,  the  tom  w£/2R  leads  exclusively  to 
tho  quantity (-1) in  the  coefficient  of  $  in  tho  a  term  of  tho  determinants! 


equation  (39).  The  quantity  persists,  and  appears  in  the  denominator  of 
the  final  simplified  expression  for  the  critical  pressure  of  a  reinforced 
cylinder,  (48). 

For  tho  caso  of  a  cylinder  with  simple  end  supports,  which  amounts 
to  setting  tho  ring  parameter  torra  of(4 8)  equsl  to  aero,  tho  simplified 
equation  of  von  kiisos  is  obtained  with  tho  *~i  retained  in  the  denominator, 
von  Mises  then  drops  the  *-l  ih  comparison  with  terms  m2  and\zin  order  to 
compensate  partially  for  his  previous  approximations* 

Thus,  if  the  term  w2/2R  is  not  used,  the  effect  would  be  that  the  ~1 
in  tho  denominator  of  (48)  would  novar  appear,  so  giving  (on  suitable  modi¬ 
fication)  tho  final  simplified  von  Mises  result  directly.  Nevertheless,  as 
is  woll  known  tho  approximate  von  Mises  result  does  not  roduco  to  tho  Bryan 
solution. 

It  is  of  intorest  to  noto  that  tho  value  of  the  critical  pressure  as 
obtained  from  the  Donnell  equation  for  circular  cylindrical  ahells  agrees 
with  the  simplified  von  Misos  result;  the  -1  town  does  not  appear  in  tho 
denominator. 
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APPENDIX  I 


Discussion  of  Approximations  ms  do  In  Pibal  number  107. 

1.  The  complete  expansion  of  the  general  determinant  (39)  leads  to  a 

cubic  equation  in  $  ,  $  and  k.  As  an  approximation  the  equation  was 
1  2 

linearized  and  took  the  form 

C+Cks*n<5+C$  . 

J  2  3  ,  4g 

Wo  now  discuss  the  validity  of  neglecting  the  terms  containing  powers  of 

$  ,  $  and  k  of  order  higher  than  unity. 

1  2 

By  definition 

k  »  (h/R )2/l2  j 

and 

<2  -2<5  -  Rp/K  «  (p/E)(H/h)U~V8)  . 

i  2 

By  definition  of  a  thin  shell,  tfa  thicKnose  h  oust  bs  small  compared 
with  the  radius  R,  i.e.  (h/R)  «  1.  Honco 

k2  -  (h/H)yi44  «  k  «  1. 

Similarly 

k  3  «  k  2 

The  buckling  pressure  p  must  always  bo  small  compared  with  the 
module!  of  elasticity  E  of  the  material  of  which  the  shell  is  constructed. 
Also,  for  the  range  of  shell  dimensions  considered  in  this  work,  the 


p 


38 


ratio  (p/r)(R/h)  «  1*  In  tho  numerical  uxaj-iplu  quoted  in  Pibal  number  167 
the  f  i-'rurc-s  were 

h/R  «  0, 006068  , 
pcrit  =  y7°  p,5,i  * 

together  with 

k  =.  3*068  x  10“°  , 

_3  _ 

$  =3  4*847  x  10  ,  (from  48), 

1 

^  =  2  •  424  x  1?  ^  , 


which  Justifies  the  linearisation  for  this  caso.  similar  results  follow 
for  the  ranfos  of  shell  and  frame  geometry  which  formed  the  bases 
of  oharts  presented  In  subsequent  reports. 

2.  In  the  bony  of  the  report  it  was  statud  that  tho  t-nns  w  ux/2R, 
v2/4  andu2/4  in  (2o)  when  integrated  lee  to  torus  of  order  hi-iier  than 

X  X 

those  given  by  the  ro:.»indor  of  tho  terms  under  the  integral*  Consequently 
they  ware  ignored. 


Consider  the  t-rm  w  uv/2R  in  (20).  It  can  be  shown  that  this  quantity- 

leads  only  to  the  turn  $  \  /2  in  the  cioa-ents  a  arid  a  of  tho  genrrsl 
_  1  13  31 

determinant  (39).  Upon  expansion  of  the  determinant  this  tom  contributes 

to  C  .  that  is  the  coefficient  of  $  .  The  minors  of  the  terms  a  and  a 
3  ’  1  13  31 

show  that  the  contributions  to  by  this  term  are  of  order  m\  ,  m4\k, 
m2\3  k,  with  ether  quantities  of  lower  order  in  m  and  \  .  On  the  other 
hand,  uxa.-dnation  of  the  whole  expression  for  C3  ,  (44c)  of  Pibal  number  167, 
8 how s  that  turns  exist  of  order  m2  X  ^ ,  m1*-  \2  y  ir.6  , 


y 


with  othor  tents 
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ratio  (p/c)(R/h)  «  1*  In  t’.io  numerical  example  quoted  in  Pibal  number  167 
the  fi'rures  were 

h/R  =  0,006068  , 

P„J^^  a  970  p.S.i  , 

together  with 

k  =  3*068  x  10“°  , 

-3  _ 

$  a  4*847  x  10  ,  (from  43), 

1 

$  =  2*/,24  x  lV*^  , 

2 


which  justifies  the  linearisation  for  this  csbo.  bimilar  results  follow 
for  the  ranges  of  shell  and  frarao  geometry  which  formed  the  bases 
of  charts  presented  in  subsequent  reports. 

2.  In  the  body  of  the  report  it  was  statud  that  the  t-rms  w 

v2/4  and  u2/4  in  (20)  when  integrated  lec  to  terus  of  order  hiher  than 
X  X 

those  given  by  the  remainder  of  tho  terms  under  the  integral.  Consequently 


they  ware  ignored. 

Consider  the  t*.rm  w  uv/2R  in  (20).  It  can  be  shown  that  this  quantity 

leads  only  to  the  term  $  \/2  in  the  elements  a  and  a  of  tho  general 

1  13  31 

determinant  (39).  Upon  expansion  of  tho  determinant  this  tom  contributes 

to  C  .  that  is  the  coefficient  of  $  .  The  minors  of  tho  terms  a  and  a 
3  *  ]  13  31 

show  that  the  contributions  to  C  by  this  term  are  of  order  ,  m4  \  k, 

3 

k>  with  ether  quantities  of  lower  order  in  m  and  \  •  On  the  other 
hand,  examination  of  the  whole  expression  for  ,  (44c)  of  Pibal  raintour  167, 
shows  that  terms  exist  of  order  m2\  **  ,  r\**  \2  ,  ni6  ,  with  other  terms 


> 


p 
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in  m  and  X  of  lower  order.  It  will  bo  shown  in  tho  following  that  the 
latter  terms,  m 2  lA  ,  m^X2  ,  m6  ,  are  alone  dominant  in  C  3  . 

Similar  considerations  will  also  show  that  the  terms  v2/ /, ,  U2/4 

J..  ^ 

lead  to  terms  in  C  of  lmwr  order  in  m  and  X  than  the  so-oallod  dominant 
3 

terms. 

To  show  that  the  higher  powers  of  m  and  X  dominate  in  the  expression 
for  ,  and  ind-.od  in  the  expressions  for  C}  ,  Cg  >  and.  as  well,  wo 
recall  the  dej’initions  of  iu  and  X  .  Thus,  m  is  tho  motor  of  cireumf orontial 
half-waves,  while  X  »  (nnR)/L,  where  R/L  is  th6  radius  of  tho  shell  divided 
by  the  distance  between  reinforcing  rings,  and  n  ia  tho  number  of  longitudinal 
half-waves. 

For  shells  of  the  typo  considered  hero,  i.e,  with  closely  spaced 
rings,  tho  A/L  ratio  is  of  the  ordor  of  2  to  (>.  In  tha  •oproasion  for  C? 

(44a)  of  Pibal  number  167)  the  first  two  terms  on  the  ritfit  hand  side  con¬ 
tain  X4  terms,  while  the  remaining  terms  contain  \  raised  to  a  lower 
power.  Consequently,  all  but.  the  first  two  terms  can  bo  neglected. 

Similar  considerations,  and  the  foct  that  k«  1,  indicate  that  for 
C2  ,  C^and  C, ,  all  terms  other  than  tho  first  can  be  neglected. 

For  sholls  with  large  R/L  ratio  it  would  be  reasonable  to  expect 
the  number  m  of  circumferential  lobes  to  be  lar/re.  Thus  terms  containing 
tha  higher  powers  of  >n  would  dominate  in  the  expressions  for  tho  coeffi¬ 
cients  C,  ,....,  C ,  .  It  so  happens  that  tho  terms  already  neglected  in  vir- 

tuo  of  the  foregoing  discussion  are  also  the  t..nos  of  lower  pemure  of  m  than 
those  r;tainod. 

V!e  conclude,  therefore,  that  the  final  expressions  for  C  C 

1  4 

which  can  safely  be  used  for  computing  the  critical  proseuie  are  as  given 
by  (47)  of  Pibal  number  167.  'tote  that  these  expressions  sro  only  valid 
for  larg?  vo Ivies  ol'  m. 


p 
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APPjSKDIJC  II. 


Dlacujsion  of  Approximations  made  to  the  Kotinard  expressions. 

1.  In  the  body  of  the  report  it  was  stated  that,  correspond*  ng  to  neglect 

of  the  Pibal  terns  w  u^/R,  v^/4,  u^/4,  the  Kennard  terras  u  w^/R  and 

v  7i  /R  can  likewise  be  negloeted. 

6 

Consider  the  tern  u  w^/R  in  the  integral  of  (19).  As  was  stated 

earlier,  this  term  lends  to  quantities  in  the  coefficient  C  ?  which  are 

negligible  compared  with  the  results  of  other  tdrms  in  that  integral. 

It  ' in  easily  be  seen  that  the  quantity  u  w  /R  leads  to  the  term  $  A 

x  I 

in  the  elements  b  ^  ^  and  b3)  of  the  general  determinant  (24).  As  in  Appendix 

1,  we  again  consider  the  minors  of  these  two  elements.'  They  lead  to  termB 
of  order  m2,\  ,  m^Ak,  m2\3k  together  with  lower  order  terms.  Examin¬ 
ation  of  the  total  expression  for  C  ? ( see  balow)  shows  the  dominant  terms  to 
be  of  order  m2A  **  ,  m^A2  and  m6  .  to  conclude  that  omission  of  tho 
term  u  wx/R  negligibly  affects  the  buckling  pressure. 

2.  It  was  shown  in  the  body  of  tho  report  that  tho  quantity  3>(  in  the 
elements  b^  and  of  tho  determinant  (24)  must  be  neglected  in  order  to 
arrive  at  Bryan's  result  for  tho  buckling  of  a  cylinder  of  infinite  length. 
This  quantity  can  bo  shown  to  result  from  tho  term  v  wg/R  in  the  integral  of 
(19).  However,  even  if  the  turn  was  retained  it  would  only  lead  to  non¬ 
dominant  terms  intho  expression  for  . 

li.  The  expressions  for  C5  ,  C  ^  and  G?  ,  resulting  from  the  expansion  of 
tho  determinant  (24  ),  arts* 

C  «  t(lf  v)/2]  f  A^  (1-  v2)  +  R  (m2  +  a2)2  +  R  R.(  A2+ 
p  e  O  o  i*  1— V 

+  X2R,  +  2  A  (m2  -  vA^R^-  R^t  A2*  2m/ ( 1-  v)]j  f 
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Cfi  -  C(l-v)/2]  f(m?-+\2)4  *  2(m  2  +X2)2(  1-m2  - vX2)  -4(k2  *  X2)(n2X2)(l~  v  ) 
-  (a2  +  vX  2)2+\2(l-4m2)RA  ♦  R;(X2  +  f~>COn2+X2)2  +  l-2m2  -  SvX2} 

■  +  2  X  (2m4+  2m2  X2  ~  a2  ~  vX2  )  R  -  X2R2 

*iC  AC 

+  Rct(n2  +  X2)2  +  X^]}  » 

C?  -  t(l- v)/2]  |  (m2  -I  +X2/2)(m2  +  X2)2  -  2  X2  (a2  +  vx2) 

+  k(m2 -1  +X2/2)(mz+  X2)2  -  2m2  X2  (2m2  +  %  \2_i) 

♦  [2/0- V  >][»*♦  (l~v)(l  +  k)][  RA(m2-l+  ♦  2  R^lj  . 


The  various  simplifications  of  those  coefficients  ere  the  seme  as 
those  outlined  in  §  2  of  Appendix  I. 


p 
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APPBSDTX  III. 

Remaining  Pi ffc micas  bqtvtuen. Pibal  number  167  "'id  the  present  report. 

It  m3  stated  in  the.  body  of  thu  report  that  ossent iaily  the  differences 
between  Pibal  number  16?  and  the  present  report  rociuce  to  showing  that  (19) 
and  (20)  are  virtually  the  same  ntmaricaHy,  that  is 

(19)  -  (20)  -  0  . 

After  accounting  for  several  of  the  terms  in  this  equation  the  remaining 
difference  D  was  shown  to  be 


D  -  -  (l/2)Rp  J‘  (wf  +  w  w  -  2r/R)  dA  , 
8  88 


W  l  , 

*  -  (l/2)Rp  J  f  (w*  +  w  »»88  ~  2w/K)  dx  ds.  (a) 


The  deflection  function  w  assumed  in  Pibal  number  167  was 


w  »  C  cos  £L_  sinAS  .  (b) 

R  L 


Consider  the  torn  ~  on  the  right  hand  side  of  (a).  Using  (b)  it  can  bo 
R 

rewritten  as 

2jtR  L  2nR  L 

f  S  ("/&)  dx  ds  -  |  f  /cos  SSL  sin  (  Xx/l,)dx  ds  •  0  , 

o  o  «  o  o  R 


This  term  would  also  b»  zero  if  w  woru  any  deflection'  function  of  the 

form 


or 


w  ■  X(x)  cos  (ms/R) 
w  »  X(x)  sin  (ms/R) 


(c) 


where  X(x)  is  any  suitably  continuous  function  of  x 


p 


p 
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next,  consicJ er  the  remaining  two  terms  of  (a).  For  the  asaunwd 
function  (b)  wo  havo 

£  w^  dA  -  t-  J*  v*  i»  dA 
7»  8  A  88 

so  that  tha  two  terras  cancel  each  other.  This  result  is  also  true  for 
any  other  deflection  functions  w  of  the  form  (c). 

We  conclude,  therefore,  that  for  any  w  of  the  form  (c),  the  final 
results  of  this  report  coincide  numerically  with  the  final  results  given  in 
Pibal  nuutoer  167* 
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.APPamiz  iv 

Buck  I:;  n  g  of  a  Cylinder  tinder  Uniform  Pressure 
Note  by  E.  H.  Konnord 


It  Fundamentals 

Coordinates  on  median  surface t 
x,  s  •»  a<p  ,  a  «  radius, 

<p  »  angle  about  axis. 

By  assumptions  of  Salerno  and  Levine  (PIBAL  in  equilibrium 

3tate  just  before  buckling  the  stresses  ars,  P  denoting  "axial  pressure," 

ax0  «  -  P,  aso  »  “2ap/2h  «  -  ap/h. 

In  comparison  with  those,  tho  radial  stress  is  neglected. 

Hence  equilibrium  strains  juBt  before  buckling  are 

e  «  (l^)(r-P  +  v  ap/h),  f  »  (l/E)(vP  -  ap/h), 

xo  30 

Row  the  elastic  energy  wi.il  contain  certain  terms  involving  only  the 
equilibrium  strains,  others  involving  only  the  added  buckling  strains,  and 
finally  cross  products  between  these  two  partial  strains.  Buckling  is  con¬ 
trolled  by  turns  in  tho  energy  and  in  the  external  work  (or  in  tho  "potential") 
that  ore  Quadratic  In  the  buckling  displace!.-, outs  or  their  derivatives.  Hence, 
in  general,  turns  in  the  buckling  strains  that  are  quadratic  in  this  sense 
must  bo  retained,  sineo  in  tho  cross  products  they  make  a  quadratic  coi*- 
tribution  to  tho  elastic  energy.  However,  quadratic  strains  containing 

u  (i.o.  du/d  x),  v  (  dv/dx)  or  w  need  not  be  retained,  for  squares  of 
X  s 
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these  quanfciti&s  occur  in  energy  terms  arising  fro;i,  tho  buckling  strains 

alone,  and  in  comparison  with  these  terms  any  cross  products  quadratic 

in  ux,  vgf  or  w  will  contain  also  an  equilibrium  strain  as  a  .very  small' 

factor,  so  that  such  cross  products  ore  negligible  for  practical  purposes. 

It  is  otherwise  with  rr  and  w  . 

x  s 

An  adequate  calculation  of  the  strains  after  buckling  is,  therefore, 

1/2 

e.  »  t(l+  e  ♦  u  )z  dx2  +  w2dx2]  /dx  -  1  , 

XO  A  A 

e  *t(l+B  +v  -w/a)2d82  +  ds2  ]^/ds  ~  1  , 
o  so  s  8  '  ' 


whence 


e  +  u  +  wz  /2, 
xo  x  x  '  1 


e  »  e  +  v  ~  w/a  +  w  zj2 
s  so  s  '  & 


Let  U  -  extensional  elastic  energy,  E  *  E/(l  ~  v2). 

*  8 

Then 

ui  "  ^3 /2)JV(e,8+  e2  +  2v  c  e  )  dx  da  . 

1  3  x  s  x  s 

The  value  of  just  before  buckling  is 

(hB8/2)JV(ejt2+  e20  +  2V  Exq  dx  ds. 

The  increase  in  U,  by  buckling  or  A  is,  therefore,  through  terras  of  2d 
order  in  u  and  w, 

AU,  »  (hE8/2)X'>Ccxo(2ujt  +  2v  vg-2  v  w/n  +  ’* 2  +  v  w2) 

+  e  (2  vu  +  2v  -r  2v*/a  +  vw2  +  t?) 
so  x  s  x  s 

+  u 2  +  vs  +  ~  55  V/a  +  2  v  ux  vs  2  v  u^/a]  dx  ds. 


p 
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II.  Special  casu;  w  ■  0  at  ends, 


Let  .it  be  assumed  that  the  ends  of  the  cylinder  remain  piano. 

Then 

JX  ux  dxds  mj*  (u^  »-  uQ)ds,  Uj.  »  u  at  x  «  L  , 

u0  »  u  at  x  1  0  , 


Also, 

So  v  dxds 
s 


2  na 

J  dx  f  v^da 


0  . 


Hence,  substituting  for  e  and  e  , 

'  XO  XB  » 


AU?  «  (hEs/2)XT(u2  +  v|  +  *2/a2  +  2  v«x  -  2vuy  w/a  | 

*-  2  V8?i/a)dxds  +  ap,r(«/a  ~  ti|/2)dxde 

zt& 

-  ~  wX"  Vi ^  dxds  -  I'lP  s  -  uQ)da. 


(Of) 


[For  comparison  ?.:Lth  PIBAL  167,  change  notation  thus: 


.a  E  , 


[Term  in  E  equals  1st  3  terms  of  PIBAL  l67-(3);  last  term  in  (3)  is 

shear  energy,  not  hero  considered.  Terms  in  p  and  P  in  Auj  represent 

cross-product  terms  vthich  are  on&tted  in  PIBAL  167.3 

Work  by  end  forces.  If  n  ■  0  at  the  ends  (but  not  otherwise,  since 

2d~order  terms  in  the  displacement  must  be  included) 

ana  ,  , 

W,  -  -  hP  S  (ur  **  u  )  ds  {r' 

*  0  "  0  ' 

(providod  P  is  uniform  around  the  shell). 
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Here  P  .is  assumed  constant  during  the  displacement. 

[P2BAL  167  hsa,  from  (20), 

W  -  -  V(  «  OiP/2)//  (ux2  +  v  2  +  wjjdxda.  (/S') 

The  deduction  of  (20)  is  not  understood.] 

Work  by  radial  pressure  p. 

This  work  equals  p  times  decrease  in  volume. 

If  w  *  0  at  the  ends,  effects  of  end  dieplacotacnts  can  be  ignored 
as  arising  from  thirdr-order  changes  of  volume,  since  FG  and  F’G*  are  them¬ 
selves  of  2d  order  in  u,  w,  (Fig.  j)  • 

Any  radius  is  shortened  by  [w]^  or 

w  -  u  w  -  v  w 
x  a 

Hence  decrease  in  volume  is,  through  quadratic  terms, 

2na 

l/r>7ia2dx  ~S dx  S  -5(3  -  «  +  ui*x  +  vw8)[(a  -  «  +  u  "yt-vw^j/ajda 
o  ^ 

•f  to.  j  (w  -  »a/2a  -  U  *x  -  v  *s)  ds 
and  work  is 

Wg  -  njvn  (w  -  wV2»  ”  u  \  -  v  wa)dxds  .  M 

[FIBAL  167  has  in  (25) 

W  »  PjJX*  -  w2/2a  +  u  w  -  aw  »  /2)dxds,  fo/' 

The  deduction  is  not  well  understood.] 
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Overall  comparison  with  PIBAL  167. 

For  an  overall  comparison,  assume,  as  in  tho  applications,  that  P 
is  due  to  p  on  end  closures,  so  that  (noarly  enough) 

P  -  ap/2h  . 

Write  (  AU,)  for  the  part  of  A  U,  in  (  a  )  not  containing  P  or  p,  this 
b 

part  agreeing  with  PIBAL  167. 

Then,  from  (  a  )  (  0  )  (  y)« 


A U  r.  W  -  W  -( AU  ), +ap//(-Sf  -  5L  „  OL  +  'HO*  +  !!s)dxds  (O 
1  1  *  1  »  2a2  4  2  a  a 


[PIBAL  167  would  get,  from  &U,  «  (  AU,  )  and  (  y*)* 

40,  UB,)b 

2  w 

+  (_  JL  +  __52  „  uxw/a  -  u2  -  v2  „  »f)dxds. 


The  first-order  (l)  term  appearing  hrare  is  eliminated  by  the  special 
Zn 

assumption  that  J*  wd<p  »  0  .] 


4 


p 
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FIG.  i  COORDINATE  SYSTEM  AND 
DISPLACEMENTS  FOR  CIRCULAR 
CYLINDRICAL  SHELL 


FIG.  2  AXIAL  DISPLACEMENTS  OF 
LONGITUDINAL  ELEMENT  OF  SHELL 


FIG.  3  DISPLACEMENT  OF  SHELL 
IN  LONGITUDINAL  PLANE 


